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Abstract 

This  investigation  is  concerned  with  the  diffusion  of  an  axial 
load  from  a  bar  of  arbitrary  uniform  cross-section  that  is  immersed 
in,  up  to  a  finite  depth,  and  bonded  to  a  semi-infinite  solid  of  distinct 
elastic  properties.  The  bar  is  perpendicular  to  the  plane  boundary  of 
the  embedding  medium.  The  determination  of  the  desired  resultant 
force  in  the  submerged  bar- segment  is  reduced  to  a  Fredholm  integral 
equation  by  means  of  an  approximative  scheme  developed  and  tested 
earlier  in  connection  with  a  more  elementary  three-dimensional  load- 
transfer  problem.  Extensive  numerical  results  illustrating  the  decay 
of  the  bar-force,  appropriate  to  various  choices  of  the  governing 
geometric  and  material  parameters,  are  presented  for  the  particular 
case  of  a  bar  of  circular  cross-section. 

Introduction 

Two-dimensional  load-transfer  problems,  such  as  those 
occasioned  by  the  diffusion  of  load  from  an  elastic  stringer  into  a 
£ 

The  results  communicated  in  this  paper  were  obtained  in  the  course 
of  an  investigation  conducted  under  Contract  Nonr-220(58)  with  the 
Office  of  Naval  Research  in  Washington,  D.  C. 
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coplanar  elastic  sheet,  have  attracted  repeated  attention,  chiefly 
because  of  their  relevance  to  the  stress  analysis  and  design  of  air¬ 
craft  structures.  The  rather  extensive  literature  on  problems  in  this 
category  goes  back  in  its  origins  to  a  fundamental  paper  by  Melan  [l], 
in  which  he  dealt  with  the  transfer  of  an  axial  load  from  an  infinite 
stringer  to  an  all-around  infinite  elastic  sheet*.  The  problem  con¬ 
cerning  the  transmission  of  an  axial  load  from  a  transverse  stringer 
(tension-bar),  a  finite  segment  of  which  overlaps  with,  and  is  contin¬ 
uously  bonded  to,  a  semi-infinite  elastic  sheet,  was  initially  posed  by 
Reissner  [2].  References  to  related  publications  are  to  be  found  in 
[3],  [4],  [5],  which  contain  some  recent  contributions  bearing  on 
Melan's  and  Reissner's  problems. 

The  diffusion  of  load  from  a  bar  into  a  three-dimensional 
elastic  medium  presents  appreciably  greater  analytical  difficulties. 

A  comparatively  simple,  if  somewhat  artificial,  example  of  this  kind 
is  supplied  by  the  spatial  analogue  of  Melan's  problem,  investigated 
in  [6],  There  we  considered  an  infinite  cylindrical  bar  fully  im¬ 
mersed  in,  and  continuously  bonded  to,  an  elastic  medium  of  distinct 
mechanical  properties  occupying  the  remainder  of  the  space;  the  bar 
was  assumed  to  be  subjected  to  an  axial  loading,  confined  to,  and  uni¬ 
formly  distributed  over,  one  of  its  cross-sections.  Our  objective 
consisted  in  determining  the  decay  of  the  resultant  bar-force  induced 
by  the  given  loading.  This  task  was  performed  rigorously,  within 
three-dimensional  elastostatics,  for  the  special  case  of  a  bar  of  cir¬ 
cular  cross-section.  The  results  thus  obtained  were  then  compared 

Within  Melan's  approximate  formulation,  this  problem  is  mathe¬ 
matically  equivalent  to  the  one  in  which  the  stringer  is  attached  to 
the  edge  of  an  elastic  sheet  occupying  a  half-plane. 
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and  found  to  be  in  favorable  agreement,  with  those  deduced  in  [6]  by 
an  approximate  method  that  is  applicable  to  a  bar  of  arbitrary  cross- 
section. 

The  work  summarized  in  [6]  was  undertaken  primarily  as  a 
pilot  study  for  the  present  investigation  of  the  three-dimensional 
counterpart  of  Reissner's  problem,  which  may  be  described  as  follows. 
A  finite  or  semi-infinite  elastic  bar  of  arbitrary  uniform  cross-section 
is  immersed  in,  up  to  a  finite  depth,  and  continuously  bonded  to,  an 
elastic  half-space  of  different  material  properties,  the  cylindrical 
boundary  of  the  bar  being  perpendicular  to  the  plane  boundary  of  the 
half- space.  The  applied  loading  is  confined  to  the  projecting  portion 
of  the  bar  and  is  taken  to  be  statically  equivalent  to  a  centroidal  axial 
force.  We  seek  the  decay  of  the  resultant  bar-force  in  the  embedded 
portion  of  the  bar,  as  a  measure  of  the  diffusion  of  the  applied  load 
into  the  surrounding  semi- infinite  medium.  This  problem  is  one  of 
basic  importance  in  structural  engineering  and  soil  mechanics  because 
of  its  interest  in  connection  with  anchor  bars  and  pile-supported  build¬ 
ings.  Moreover,  its  solution  is  not  readily  accessible  by  rigorous 
means  even  for  a  bar  of  circular  cross-section. 

The  approximative  scheme  employed  by  Melan,  and  later 
adopted  by  Reissner,  in  dealing  with  the  plane  load-transfer  prob¬ 
lems  considered  in  [l]  and  [2],  respectively,  has  no  direct  analogue 
in  the  analysis  of  the  corresponding  spatial  problems.  To  make  this 
clear  we  recall  first  the  simplifying  assumptions  underlying  [l],  [2]. 
In  both  instances  the  bar  is  modeled  as  a  one-dimensional  elastic 
continuum,  whereas  the  sheet  is  treated  as  a  two-dimensional  elastic 
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medium  within  the  conventional  theory  of  generalized  plane  stress. 
Consequently,  both  Melan  and  Reissner  regard  the  bond-forces  exerted 
by  the  stringer  on  the  sheet  as  an  ideal  line-load  and  impose  as  a  bond 
condition  the  requirement  that  the  axial  strain  in  the  bar  match  the 
appropriate  extensional  strain  in  the  sheet.  This  condition,  with  the 
aid  of  the  familiar  two-dimensional  singular  solutions  pertaining  to  a 
concentrated  load  at  an  interior  point  of  an  infinite  or  semi-infinite 
elastic  sheet,  leads  to  a  singular  integro-differential  equation  for  the 
desired  bar-force,  the  integral  thus  emerging  being  convergent  in  the 
sense  of  its  Cauchy  principal  value. 

A  strictly  one-dimensional  bar-model,  in  conjunction  with  a 
bond  condition  based  on  the  assumption  of  ideal  line-contact,  is  no 
longer  feasible  in  dealing  with  the  transfer  of  load  into  a  three- 
dimensional  elastic  medium:  here,  in  contrast  to  the  situation  encoun¬ 
tered  in  two  dimensions,  the  strain  field  due  to  a  line-load  is  unbounded 
along  the  line- segment  of  load-application^.  Indeed,  the  formal  ana¬ 
logue  of  Melan's  formulation  in  the  present  circumstances  gives  rise 
to  a  bond  condition  that  involves  a  divergent  integral  and  is  therefore 
devoid  of  meaning. 

To  circumvent  the  difficulty  described  above,  it  is  essential  to 
amend  Melan's  assumptions.  A  suitable  approximate  formulation  of 
the  present  problem  is  given  in  Section  1.  This  formulation  is  based 
on  an  adaptation  of  the  approximative  scheme  devised  in  [6],  which 
was  in  turn  suggested  by  a  corresponding  treatment  of  Reissner's 


1) 


Recall  that  the  strain  singularity  associated  with  a  concentrated 
load  in  three- space  is  one  order  higher  than  its  two-dimensional 
counterpart. 
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plane  problem,  carried  out  in  [5]*.  The  remainder  of  the  section  is 

devoted  to  reducing  the  determination  of  the  desired  resultant  force  in 

the  bar  to  the  solution  of  an  integral  equation  of  Fredholm's  serond 

kind.  An  essential  tool  for  the  analysis  in  Section  1  is  supplied  by 

Mindlin's  [7]  solution  to  the  problem  of  a  concentrated  load  acting  at 

2 

an  interior  point  of  an  elastic  half- space  . 

In  Section  2  we  specialize  the  characterization  of  the  bar-force 
deduced  in  Section  1  for  the  particular,  physically  most  significant, 
case  of  a  bar  of  circular  cross-section.  Here  we  also  establish  cer¬ 
tain  asymptotic  properties  of  the  solution  appropriate  to  the  circular 
bar  and  present  illustrative  quantitative  results  based  on  the  numerical 
solution  of  the  Fredholm  integral  equation  that  governs  the  diffusion 
of  the  applied  load  in  the  present  instance. 

1.  Approximate  formulation  of  problem.  Reduction  to  a  Fredholm 
integral  equation. 

We  proceed  now  to  an  approximate  mathematical  formulation  of 
the  three-dimensional  analogue  of  Reissner's  problem  described  in  the 
Introduction.  To  this  end,  let  {Oix^x^.x^}  be  a  rectangular  cartesian 
coordinate  frame  (Figure  1)  with  the  origin  O,  spanning  the  three- 
dimensional  Euclidean  space  E.  Further,  let  x=(x^,  x^,  x^)  denote  the 
position  vector  of  points  in  E  and  call  e  the  unit  base-vector  in  the 
Xj -direction. 

Consider  a  finite  or  semi-infinite  cylindrical  elastic  bar,  the 
longitudinal  centroidal  axis  of  which  coincides  with  the  x^-axis.  We 

^  See  Section  3  of  [5]. 

2) 

■  See  Mindlin  [8l  for  a  systematic  derivation  of  the  solution  to  his 
problem. 
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call  II  the  generic  open  cross-sectional  region  of  the  bar  and  assume 
II  to  be  a  bounded  plane  region,  whose  boundary  ^  dll  is  a  simple  closed 
curve,  consisting  of  a  finite  number  of  arcs  with  continuous  curvature. 
Suppose  H  to  be  the  open  half- space  given  by 


H={x|x€E,  x,>0} , 

~  ~  J 

(1. 1) 

designate  by  D  the  (cylindrical)subdomain  of  H  defined  by 

D={x |  (Xj,x2)€II,  0<x3<t}, 

(1.2) 

and  set 

nz={x|(xr x2)€n,x3=z}  (osx^t), 

(1.3) 

where  11^  is  the  open  cross-section  of  7  located  at  x^=z.  In  particular, 

n  and  II.  stand  for  the  interiors  of  the  terminal  cross-sections  of  15, 
o  •o 

situated  at  x^=0  and  x^-l,  respectively  (Figure  1).  Next,  consider  an 
elastic  body  ("embedding  medium")  whose  interior  occupies  the  domain 
H-E  and  suppose  one  end-portion  of  the  bar  occupies  E  and  is  continu¬ 
ously  bonded  to  the  surrounding  solid  throughout  the  interface  BD-I^. 
Accordingly,  D  represents  the  interior  of  the  embedded  segment  of  the 
bar  and  l  is  the  length  of  this  submerged  segment.  Finally,  assume 
the  entire  loading  is  confined  to  the  projecting  portion  of  the  bar  and  is 
statically  equivalent  to  a  force  -epQ  acting  along  the  x^-axis,  so  that 
Po>0  corresponds  to  a  tensile  loading.  We  seek  the  resultant  axial 
bar-force  p(x^)  (Ofix^st),  which  governs  the  load-diffusion  into  the  em¬ 
bedding  medium,  on  the  assumption  that  the  latter,  as  well  as  the  bar 
itself,  are  homogeneous  and  isotropic  linearly  elastic  solids  of  possibly 
distinct  material  properties.  Let  r\  and  v,  in  this  order,  stand  for 


TT 


If  S  is  a  point-set  in  a  Euclidean  space  of  two  or  three-dimensions, 
we  write  3  and  8S  for  the  closure  and  the  boundary  of  S,  respectively. 
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Young's  modulus  and  Poisson's  ratio  of  the  embedding  medium; 

Young's  modulus  of  the  bar*  will  be  designated  by  r/ . 

The  present  approximate  treatment  of  the  problem  aims  at, 
and  is  confined  in  its  physical  relevance  to,  a  bar  with  a  maximum 
cross -sectional  diameter  that  is  suitably  small  compared  to  the 
length  of  the  submerged  bar- segment.  As  a  first  step  we  extend  the 
embedding  medium  throughout  the  half- space  H  and  consider  an  elastic 
body  B  in  17  with  the  elastic  constants  t|,  v  of  the  original  material  in 
H-D.  Next,  we  introduce  a  fictitious  reinforcement  of  the  body  B 
throughout  the  region  D  {Figure  2).  The  reinforcement  B^  is  chosen 
in  such  a  way  that  the  "composite  solid"  occupying  D  is  "equivalent" 
to  the  actual  immersed  bar- segment  B*  in  the  sense  that  its  extensional 
stiffness  is  the  same  as  that  of  B*.  We  thus  assign  to  the  reinforce¬ 
ment  B^  the  modulus  of  elasticity  r\ ^  given  by 

1VT,"T1*0-  (1.  4)2 

In  what  follows  we  treat  the  extended  embedding  medium  B  as  a 
three-dimensional  elastic  continuum  within  the  framework  of  classical 
elastostatics.  In  contrast,  we  regard  the  fictitious  reinforcement  B^ 
as  a  one-dimensional  elastic  continuum  as  far  as  its  constitutive  law 
and  equilibrium  conditions  are  concerned.  Accordingly,  B^  is  governed 
by  the  stress-strain  relation 

>fP*(z)=1V*<z>  (1.5) 

if  A  is  the  area  of  FT,  p^z)  is  the  scalar  axial  force  in  B#  at  x^=z, 
positive  if  tensile,  and  e^z)  is  the  associated  axial  strain.  Further, 

*)  Poisson's  ratio  of  the  bar  does  not  enter  into  the  subsequent  analysis. 

^  Note  that  (1.4)  rules  out  the  practically  uninteresting  case  in  which  V<r|. 
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the  equilibrium  requirement  for  B+  furnishes  the  differential  equation 

p^(z)+q(z)=0  (0<z<t),  (1.6)1 

where  £q(z)  is  the  "bond-force"  per  unit  bar-length  exerted  by  B  on  B* 
at  x^=z.  In  addition  to  the  distributed  bond-force  of  scalar  density 
q(z)  (0<z<t),  B+  is  subjected  to  the  external  axial  forces  “®p„,(0)  and 
ep^t),  concentrated  within  the  terminal  cross-sections  IIo  and 
respectively:  -fcp^O)  is  the  portion  of  the  applied  load  -£pQ  trans¬ 
mitted  to  B^  directly,  whereas  ep^-t)  is  the  bond-force  communicated 
by  B  to  B^  at  x3=i  (See  Figure  2). 

The  forces  exterral  with  respect  to  the  body  B  (extended  matrix) 
consist  of  the  distributed  bond-forces  of  lineal  density  -£q(z)  (0<z<£) 
exerted  by  B^  on  B,  together  with  the  load-portion  -g[p  -p^O)]  trans¬ 
mitted  to  B  directly  at  x3=0  and  the  bond-force  concentrated 

within  TT^.  We  assume  the  foregoing  forces  to  be  uniformly  distributed 
over  the  respective  cross-sections  IT  (0<z<t),  II  ,  and  II.  (See  Figure  2). 

Let  a. j  and  e„  (i,  j=l,2,  3)  be  the  cartesian  components  of  the 
fields  of  stress  and  strain  induced  in  the  body  B  by  the  system  of  exter¬ 
nal  forces  acting  on  B.  Further,  for  every  x€R-!I^  and  every  C€[0,-t], 
let  6  ..(g,  (!)  and  e..(x,  C)  represent  the  stresses  and  strains  at  the  point 
£  of  the  semi-infinite  medium  B  due  to  a  uniform  body-force  distribu¬ 
tion  over  the  closed  disk  11^,  acting  in  the  negative  x^ -direction,  t>  e 
resultant  applied  force  having  unit  magnitude.  With  the  aid  of  these 
influence  fields  one  has,  for  every  point  x€H-fl  -IT^, 


aij<2S)=CPo-Pe<0)]dij^ ^J^Ctfijlx. C)dC.  I  (1>7) 


If  v  is  a  differentiable  function  of  one  variable,  we  write  v  for  the 
derivative  of  v. 
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(1.7) 

cont. 


Integral  representations  for  the  influence  fields  and 

entering  (1.7)  are,  in  turn,  immediately  obtainable  from  Mindlin's  [7], 

[8]  solution  for  a  concentrated  load  applied  at  a  point  of  an  elastic 

half-space.  Indeed,  designating  by  2r..(x,  C)  and  e..(x,  C)  the  stresses 

^  J  ^  J 

and  strains  of  Mindlin's  solution  appropriate  to  the  semi-infinite  solid 
B  and  corresponding  to  the  unit  concentrated  load  in  the  negative 
Xj -direction  applied  at  (0,  0,  C),  one  has 

fy*. C>4 Jfy*-?. C)<iAr  (1'8)‘ 


for  all  (0 Explicit  representations  for  and  e^,  in 

closed  elementary  form,  will  be  cited  later  on. 

Let  p(z)  be  the  desired  axial  bar-force  at  x^=z  in  the  immersed 
bar  segment  B*.  Then,  within  the  approximation  under  consideration, 


p  obeys 


P(z)=p+(z)+J  033(x)dA  (Oszst), 


(1.9) 


where  p(0)=p(0+)  and  p(t)=p({,-).  In  order  to  render  p  fully  determinate, 
one  needs  to  adjoin  to  (1.9)  a  suitable  bond  condition  that  links  the 
deformations  of  the  fictitious  reinforcing  bar  B^  and  those  of  the 
extended  embedding  medium  B.  For  this  purpose  we  adopt  the  require¬ 
ment  that  the  axial  strain  C*(z)  in  B*  be  equal  to  the  cross-sectional 
average  over  IIZ  of  the  extensional  strain  e^^(x)  in  B.  Consequently, 

^  Here  and  in  the  sequel  the  subscript  attached  to  the  element  of  area 
indicates  the  variable  of  integration. 
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e*^z)  =  ^J  e33^)dA  (0Sz^)* 
n 


(1.10) 


with  the  understanding  that  the  integral  in  (1.  10)  for  z=0  and  z=t  is  to 
be  interpreted  in  the  sense  of  its  corresponding  limits  as  z-*0  and  zM, 
from  the  interior  of  the  interval  [0,  £].  Equations  (1.  5)  to  (1.  10)  suf¬ 
fice  to  characterize  p(z)  completely  over  the  range  0*z£l. 

With  a  view  toward  making  the  foregoing  characterization  of  the 
bar-force  p  more  explicit,  we  define  the  auxiliary  influence  functions 
0  and  e  through 


a(z,C)=^J  S33(x,C)dA  (0Sz*t,0£C*WO, 


g(z,  C)=4  J  e33(x,C)dA  (Q*zil,0*Q*l,zK), 

n_ 


)  (i.ii) 


(1.12) 


supplemented  by 

O(0,0)  =  O(0+,0),  0(t,l)=0(l-,l), 

e(0,  0)=c(0+,  0),  c(l,l)=z(l-,l). 

It  will  be  shown  subsequently  that  the  limits  involved  in  (1.  12)  exist. 

Further,  o(z,  C)  and  e(z,  Q,  regarded  as  functions  of  0  will  be  found  to 

possess  merely  finite  jump  discontinuities  at  C=z  for  each  fixed  z  in  the 

open  interval  (0,1).  The  physical  significance  of  0  and  e  as  stress  and 

strain  averages  associated  with  the  influence  fields  d. .  and  e..  is  imme- 
B  ij  ij 

diate  from  (1.  11),  (1.  12). 

We  now  substitute  for  033(x),  €33(x)  from  (1.  7)  into  (1.  9),  (1.  10) 
and  —  after  a  permissible  reversal  in  the  order  of  the  two  integrations  — 
invoke  (1.  11),  (1.12),  as  well  as  (1.  5),  (1.  6),  to  arrive  at 
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P(*)=P*(*)+A{[po-pi|t(0)]o(z,  0) 

l 

+P*W<J{z,t)-}p)(i(C)a(z,  C)dc}  (Osz^-t).  (1.13) 

0 

aR“P#(z)=^po-p*<0):1€(z*0) 

l 

+PJl)e(z,l)-jpJC)e(z,C)dC  (Osz^).  ( 1.  14) 1 

0 

Observe  that  the  influence  functions  0  and  e  are  computable,  for  a  bar 
of  given  cross-section  n,  from  (1.  11),  (1.  12)  in  conjunction  with  (1.8) 
and  Mindlin's  solution  [7].  Hence  (1.  14)  constitutes  an  integro- 
differential  equation  for  the  fictitious  bar-force  p^.  Moreover,  once 
p^z)  (O&z^l)  is  known,  the  actual  bar-force  p(z)  (0«z«t)  follows  from 
(1.  13). 

Our  next  objective  is  to  reduce  (1.  14)  to  an  ordinary  integral 
equation  of  Fredholm's  second  kind  and,  at  the  same  time,  to  trans¬ 
form  (1.  13)  into  a  more  convenient  form.  In  order  to  accomplish  this 
purpose  we  require  more  detailed  information  as  to  the  structure  of 
the  influence  functions  0  and  e,  which  are  related  through  (1.  8),  (1.  11) 
to  Mindlin's  [7]  solution.  The  latter  furnishes  for  every  and 

every  x€H  (x*eC), 

®33(i>C)=5^nv){2(1'v,vl,i'SC>-<x3’C)v2<2'SC) 

+2(l-v)v1(x+eC)-[(3-4v)x3+Clv2(x+eC)+2Cx3v3(x+eC)},  (1.  15) 

^  In  deducing  (1.  14)  we  excluded  the  degenerate  case  r^sO,  i.  e.  ,  T/=rt, 
which  requires  separate  attention.  We  shall  return  to  this  case  at 
the  end  of  Section  1. 
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e33<^*  C>=^7Trbjri{2(1-2^)vi<2£-*C)-(x3-C)v2<x-eC) 
+2(l-2vfvj(ic+eC)-[(3-4v)x3+{l-4v)C]v2(2+eC)+2Cx3v3(x+cC)}(  (1.  16) 


where 


8x”  Mil 


We  now  define  functions  U  and  Vn  (n=l,  2,  3)  by  means  of 


U<5>=-*J 


FH 


for  all  x€E, 


(1.  17) 


(1.  18) 


Vn(2,=2“nU(2)  for  an  £€E-TT0  (n=l,  2,  3). 

8x3 

Thus  U  is  the  Newtonian  potential  of  a  uniform  mass  distribution  of 
density  -1/A  over  the  disk  nQ.  From  (1.  18)  follows,  for  every  ^€E-TT0, 

x3  p  dAi  if  hs-!l2-3x3  > 

vi(*)=xJ  T-T7T’  WtJ  -T--:.r"dAg, 


n 


3x.  r  5x2-3|x-?|2 

<&-r 

n  K-S  ~ 


(1.  19) 


Finally,  we  introduce  functions  Wr  (n=l,2,3)  by  setting 

Wn(z)  =  iJ  Vnte)dA  <0<lzl<oo)'  (n=l,  2,  3).  (1.20) 


As  is  readily  confirmed  with  the  aid  of  (1.8),  (1.  11)  and  (1.  15),  (1.  16), 
(1.  17),  by  recourse  to  (1.  18),  (1.20),  the  functions  a  and  e  admit  the 
representations 
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0(Zt  0='g^TT7){2(l-v)W1{z-C)-(z-C)W2(z-C)+2(l-v)W1(z+C) 
-C(3-4v)z+C3W2(z+C)+2CzW3(z+C)}  (0*zU,0*Q*t,zH),  (1.21) 

€(z*  ^)=i^T^yn{2<1-2^wi(z-C)-<*-C)w2(z-C)+2(i-2vjaw1(z+C) 

-[(3-4v)z+(1-4v)C]W2(z+C)+2CzW3(z+C)}  (O^zSt,  O^C^,  z*C).  (1.22) 

The  behavior  of  the  functions  0  and  e  is  accordingly  governed  by, 
and  may  be  inferred  from,  the  properties  of  (n=l,  2,  3)  given  by  (1.  19), 
(1.20).  The  functions  Wn,  in  turn,  are  evidently  continuous  and  possess 
continuous  partial  derivatives  of  all  orders  at  all  points  of  the  real  axis 
with  the  exception  of  the  origin.  Moreover,  (1.19),  (1.  20)  permit  one 
to  show  that 

W1(0+)  =  -W1(0-)  =  ^f,  lim[zn_1Wn(z)]=0  (n=2,  3).  (1.23)1 

z-*0 

Now  let  f  and  g  be  functions  defined  by 

' 

0( z,  C)=^[h(z-C)+f(z,  C)]  (OSz^t,  0*1*1,  zK), 

>  (1.24) 

e(z,  C)=  1 1*^  [(l-2v)h(z- CHg(z,C)  1  (0*z*l,  o*:*i,  zK), 

where 

h(z)=l  (0<z<oo),  h(z)  =  -l  (-oo<z<0).  (1.25) 

It  is  clear  from  (1.  21),  (1.  22),  (1.  23)  and  the  regularity  properties  of 
WR  (n=l,2,3)  described  above  that  f  and  g  so  defined  are  continuously 
differentiable  on  the  slit  square 

S={(z,  C)| 0*z*l,  0*Q*t,  z^C)  (1.26) 


TT 


See  the  Appendix  at  the  end  of  the  paper  for  a  sketch  of  a  proof  of  (1.23). 
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and  coincide  on  S  with  functions  continuous  on  the  closed  square 
Consequently,  the  partial  derivatives  9f/8C  and  9g/0C,  for  each  fixed 
z€[0,  l),  are  integrable  functions  of  Q  on  [0,  -t]. 

At  this  stage  we  substitute  for  0  and  e  from  (1.24)  into  (1.  13), 
(1.  14),  bearing  in  mind  (1.  12).  We  then  apply  an  integration  by  parts 
to  the  integrals  appearing  (1.  13),  (1.  14),  taking  proper  account  of  the 
jump  discontinuities  *  inherent  in  the  step-function  h.  This  procedure, 
which  derives  its  legitimacy  from  the  previously  established  regularity 
properties  of  f  and  g  —  provided  the  unknown  function  p*  (fictitious  bar- 
force)  is  assumed  to  be  continuously  differentiable  on  (0,  l)  and  con¬ 
tinuous  on  the  closed  interval  [0,  l~\  —  yields  after  an  elementary  com¬ 
putation. 


I 


(Np0.  _ 

=  2(l-v)°  ^ l-2v+g(z,  °)^  (0*z*t),  (1.27) 

l 

p(z)=^[  l+f(z,  0)]  +£  J p#(C)|7f(z.  C)dC  (0£z^).  (1.  28) 

0 

For  a  specified  cross-section  II  of  the  original  bar,  explicit 
integral  representations  of  the  functions  f  and  g  entering  (1.27),  (1.  28) 
are  obtainable  from  their  definitions  in  (1.  24)  together  with  (1.  19)  to 
(1.  22).  Equation  (1.  27)  is  therefore  an  integral  equation  of  Fredholm's 
second  kind  for  the  fictitious  bar-force  p*.  This  integral  equation  is 
readily  amenable  to  a  numerical  solution.  Once  p*  has  been  found,  the 

IT 


See  the  definition  (1.  25)  of  h.  The  range  of  integration  Co,  -tl  needs 
to  be  decomposed  into  the  subintervals  [0,  z]  and  [z,  t"]. 
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actual  bar-force  p  —  which  is  the  primary  unknown  in  the  problem 
under  consideration  —  is  determined  by  (1.  28). 

The  value  of  p(0)  may  be  extracted  directly  from  (1.  28).  For 
this  purpose  observe  from  (1.  11),  (1.  8)  that 

o( 0,  C)=0  (0<C*t).  (1.29) 

Now  (1.  29),  (1.  24),  in  view  of  the  regularity  properties  of  f  established 
earlier,  lead  to 

£cf(z*  C)lz=0=-3£f(0*  C)=0  <0<^)*  f(°+*  °>=1-  0-  30) 

Consequently,  (1.  28)  implies 

p(0)=po,  (1.31) 

so  that  —  in  contrast  to  the  situation  encountered  in  the  analogous 
treatment  of  Reissner's  plane  problem1  —  no  concentrated  load-transfer 
occurs  at  z=0. 

Equation  (1.  27)  is  inapplicable^  for  r4f{=T1  -Tl=0,  i.  e.  ,  when  the 
bar  and  the  surrounding  solid  have  the  same  modulus  of  elasticity.  In 
this  degenerate  instance  p+=0  on  [0,-t]  according  to  (1.5),  and  thus  (1.  28), 
because  of  the  last  of  (1.  30),  gives  way  to 

p 

p(z)^[l+f(z,0)3  (0<zsl),  p(0)=po.  (1.32) 

Finally,  we  note  that  elimination  of  q  from  the  first  of  (1.  7)  with 
tne  aid  of  (1.  6),  and  a  subsequent  integration  by  parts,  lead  to  the 
representation 

l 

V~)=poV~  0)+Ip*(C,l?dij(~’ C)dC  for  al1  (1- 33) 


^  See  Section  3  of  [5]. 

2) 

"  See  the  footnote  attached  to  (1.  14). 
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of  the  stress  field  in  the  original  embedding  medium,  the  influence 
field  being  known  from  (1.8)  and  Mindlin's  solution  [7]. 


2.  Specialization  to  a  bar  of  circular  cross  -section.  Asymptotic  and 
numerical  results.  Discussion." 

In  this  section  we  specialize  the  characterization  (1.  27),  (1.  28) 
of  the  bar-force  p  by  considering  the  particular  case  of  a  bar  of  cir¬ 
cular  cross-section.  Thus,  henceforth 

Il  =  {(x1,x2)|xj+x2<a2),  A=ira2,  (2.1) 

where  "a"  is  the  radius  of  the  cross-section.  With  a  view  toward 
expressing  the  functions  f  and  g  entering  (1.  27),  (1.  28)  in  convenient 
explicit  form,  we  set 

x.=rcos0,  x,=rsin0  (0<r<oo,  Oa£0<2ir),  f 

(2.2) 

?j=pcoscp,  52  =  psincp  (0<p<oo,  0scp<2ir),  j 

and  recall  that  the  Newtonian  potential  at  (x^.x^.x^)  of  a  particle 
situated  at  (?j,  ?2*0)  admits  the  Laplace- integral  ^  representation 

C(xl-?l)2+(x2-?2)2+x3rl/2= 

oo  _ _ _ 

Jexp(-|x3|s)Jo(syi,^+p2-2rpcos(0-cp),)ds  (2.3) 

0 

for  all  &€E  except  x=(5j,  ?2»0).  Here  JQ  is  the  zero-order  Bessel 

function  of  the  first  kind.  Consequently,  the  disk-potential  U  defined 

by  the  first  of  (1.  18)  may  in  the  present  instance  be  written  as 

a  2tt  oo  _ | 

U(jc)=-— J  Jexp(-|x3|s)Jo(s/,2+p2-2rpcos(0-cp)  )pdsdcpdp  (2.4) 
Tra  0  0  0 


TT 


Watson  [9],  p.  384. 
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for  every  x€E.  Next,  recall  the  series  expansion* 

n — z - 1 

Jo(s/r  +p  -2rpcos(0-cp)  )=Jo(8r)JQ(sp) 

oo 

+2Y  J  (sr)J  (sp)cos[n(6-cp)],  (2.5) 

Li  n  n 

n=l 


in  which  J  denotes  the  Bessel  function  of  the  first  kind  of  order  n. 
n 

Substituting  from  (2.  5)  into  (2.  4)  and  bearing  in  mind  that 

a 

Jjo(sp)pdp  =  -|j1(8a)  (0<s<oo),  (2.6)2 

0 


one  obtains,  after  a  permissible  interchange  of  the  order  of  the  integra¬ 
tions  in  (2.  4)  and  upon  a  legitimate  termwise  integration  with  respect 
to  cp, 

00 

U(x)=  J  7*xp(-|jc3|8)J0(8r)J1(sa)ds  for  all  jc€E.  (2.  7)3 

0 


Equation  (2.  7)  in  conjunction  with  the  second  of  (1.  18)  now 


yields 


00 


Vn(x)=  -|[-h(x3)]n  J  sn‘1exp(-|x3|s)Jo(sr)J1(sa)ds 


for  all  x€E-8H  (n-1,2,3),  (2.8) 

where  h  is  the  step-function  introducing  in  (1.  25).  Inserting  VR  from 
(2.  8)  into  (1.  20)  and  carrying  out  the  area-integration  (over  the  circu¬ 
lar  region  II  )  under  the  improper  integral,  as  is  justified,  one  is  led 
to 


^  Watson  [9],  p.  358. 

Watson  [9],  p.  18. 

3) 

'  This  integral  representation  for  the  Newtonian  potential  of  a  uniform 
circular  disk  is  originally  due  to  Weber  [lO], 
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Wn(z)=’“^T[‘h<z):)nRn(^)  (0<lzl<oo)»  (n=l,2,3),  (2.9) 

a 

provided  R^,  for  every  positive  integer  n,  is  the  function  defined  by 

oo 

Rn(t)=an  *  J  sn  2exp(-ts)[  Jj(as)]2ds  (0<t<oo).  (2.  10) 

n  ft 


Evidently,  R^  is  continuous  and  has  continuous  partial  derivatives  of  all 
orders  on  (0,  oo). 

It  is  clear  from  (2.  9),  (1.  21),  (1.  22),  that  the  influence  functions 
o  and  e  may  now  be  expressed  in  terms  of  the  integrals  Rn  (n=l,2,3), 
the  step-function  h,  and  elementary  functions.  Once  this  has  been 
accomplished,  the  respective  "continuous  parts"  f  and  g  of  O  and  e  are 
obtainable  directly  from  (1.  24).  In  this  manner  one  arrives  at  the  fol¬ 
lowing  results  for  a  bar  of  circular  cross-section: 


f(z,C)  =  T~{(l-v)h(z-C)C2R1(|z-C!)-l3+2(l-v)R1(z+C) 

+  ^R2(lz'Cl)+[l+(3'4v)f]R2(z+C)+^TR3(z+C)5, 

a 


\ 


M2-  U) 


g(z,C)  =  (l-2v)h(z-C)C2R1(|z-C|)-l]+2(l-2v)2R1(z+C) 

+  ^R2(|z-C|)+C(l-4v)|+(3-4v)f]R2(z+C)+^R3(z+C) 

a  / 

for  all  (z,  C)  in  the  region  S  given  by  (1.  26)*.  The  required  partial 
derivatives  of  f(z,  C)  and  g(z,  Q  with  respect  to  C  are  readily  computed 
from  (2.  11)  with  the  aid  of  the  recursion  relation 


Rn<*>=4Rn+l<,>  <0<,<0o>’ 


(2.  12) 


As  was  pointed  out  in  Section  1,  the  functions  f  and  g  may  be  con¬ 
tinuously  extended  onto  S.  This  claim  may  be  confirmed  in  the 
special  case  at  hand  by  means  of  (2.  11)  and  asymptotic  estimates  of 
Rn(t)  as  t-»0,  which  will  be  cited  later  on  (See  (2.  18) ). 
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>(2.13) 


which  is  implied  by  (2.  10).  This  computation  furnishes 

|^f(z,C)  =  ^~p{(l-2v)a2[R2(|z-c!)-R2(z+C)]+a|z-C|R3(!z-C|)  ^ 

-a[C+(l-4v)z]R3(z+C)-2zCR4(z+C))  (Oszst,  OzQzt,  z^C). 

|^g(*.C)=-3{(1-4v)a2R2(|z-C|)-{1-4NH-8v2>a2R2(*+C)+aU-C|R3(U-C|> 

-(l-4v)a(z+C)R3(z+C)-2zCR4(z+C)}  (Oszst,  O^C^.  z^C). 

We  adjoin  here  the  analogous  representation  for  the  partial  derivative 
of  f(z,  C)  with  respect  to  z,  which  will  also  be  needed  subsequently. 

|jf(z,  =  {(l-2v)a2CR2(|z-C|)-R2(z+C)3+a|z-C|R3(|z-d)  ] 


(2.  14) 


-a[C-(3-4v)z]R3(z+C)+2zCR4(z+C)}  (0«z«t,  z^C). 

It  is  essential  for  numerical  purposes,  as  well  as  for  the 
asymptotic  study  of  the  actual  bar-force  p  near  the  ends  of  the  embedded 
bar- segment,  to  examine  in  detail  the  singular  behavior  at  z=C  of  the 
partial  derivatives  listed  in  (2.  13),  (2.  14).  This  behavior,  in  turn,  is 
evidently  governed  by  the  asymptotic  character  of  Rn(t)  as  t-»0.  The 


functions  Rr,  introduced  through  (2.  10),  are  expressible  in  terms  of 
the  complete  elliptic  integrals  of  the  first  and  second  kind.  Indeed,  with 
tr/2  ir/2  \ 

K(k)=  J- 


dcp 


— ,  E(k)=  j /l-k2sin2cp  dcp, 
0  /l-k2sin2cp  0 


k=k(t)  = 


2a 


>  (2.  15)1 


(0<t<CD), 


Our  previous  use  of  "E"  for  the  entire  space  ought  not  to  cause 
confusion. 
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one  has 


1 


Ri(t)4-ikCK(k)_E(k)]  (°<t<0°)» 

R2(t)  =  ^[(2-k2)K(k)-2E<k)]  (0<t<oo), 

2 

R3(t)=  "  jjjf5TK<k>+  ~  fftW'  * E(k)  (°<t<0°)» 

3 

R4(t)=feK(k)+2¥[2(a/t)2_k2]E(k)  (°<t<0°)- 


\ 


>  (2.16) 


) 


Equations  (2.  16)  confirm  the  analytic  character  of  Rr  (n=l,2,3,4)  on 

2 

(0,oo).  Now,  from  the  last  of  (2.  15)  and  familiar  asymptotic  estimates 
of  K(k),  E(k)  as  k-1,  follows 

K(k(t))=-log(t/8a)+0(t2logt)  as  t-0,  1 

>  (2.17) 

E(k(t))=l+0(t2logt)  as  t-0.  J 

Combining  (2.  16)  and  (2.  17),  one  thus  arrives  at  the  estimates 


Rl(t)  =  j+0(tl°gt)  as  t-*0, 

R_(t)=-  — [log(t/8a)+2]+0(t2logt)  as  t-0, 
R3(t)  =  £+0(tlogt)  as  t-0. 


\ 

>  (2.18) 


R4(t)=-2T+0(Iogt)  as  t-0  7 

irt 

Equations  (2.  18)  enable  one  to  remove  from  the  partial  derivatives 
(2.  13),  (2.  14),  in  closed  elementary  form,  those  contributions  that 
become  unbounded  as  C-z  (Oszst).  We  now  cite  the  decompositions  thus 
established. 

^  See  Eason,  Noble,  and  Sneddon  [ll]  for  the  first  three  of  (2.  16);  the 
last  of  (2.  16)  is  easily  deduced  from  the  preceding  one  with  the 
aid  of  (2.  12),  (2.  15). 

2) 

'  See,  for  example,  Oberhettinger  and  Magnus  [12],  p.  3. 
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C)  =  -irrr^jyir{(1-2v)logC|Z-C| /(«+C)3+irF(*.  0), 

^£g(z,  C)  =  -~((1-4v)l°g[|z-C|  /8a]-(l-4v+8v2)log[(z+C)/8a] 

+irG(z,C)}  »  z^C)» 


\ 


>  (2.  19) 


/ 


llf(z*C)=7(T^t(1_2v)lo8Clz-^/<z+c)]+wH<z*c)^ 

(0£z£{,,  OsC^-t,  z/Q, 


(2.20) 


where  F,  G,  H,  are  new  functions  to  be  discussed  presently.  Let  s'  be 
the  plane  region  obtained  by  deleting  the  origin  from  the  closed  square 
S,  so  that 


S'={(z,  C)|o<zst,0<Cstl.  (2.21) 

The  functions  F,  G,  and  H  are  bounded  and  continuous  on  S#;  moreover, 
for  any  fixed  direction  of  approach,  they  possess  limits  as  (z,  C)-*(0,0). 
Explicit  representations  are  needed  merely  for  F  and  G: 


F(z,  C)=-(l-2v){R2(|z-C|)-R2(z+C)+^log[|z-C|  /( z+C)]} 

--y{a|z-C|R3(|z^C|)-a[C  +  (l-4v)z]R3(z+C)-2zCR4(z+C)} 
a 

(O^z^t,  OzQzl,  z*C).  >  (2.22) 

F(z,  z)=i-^  +  (l-2v)[R2(2z)+^log(z/4a)] 

+  -^[(l-2v)azR3(2z)+z2R4(2z)l  (0<zs/.),  F(0,  0+)=0;  / 


G(z,C)=-(l-4v){R2(|z-d)+^log[|z-C|/8a]} 

+  (l-4\H-8v2){R2(z+Q+^log[(z+C)/8a]}--j{a|z-C|R3(|z-C|) 

,«  a 


)  (2.23) 


TT 


Recall  (1.  26). 
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-(l-4v)a(z+C)R3(z+C)-2zCR4(z+C)}  (0*zzl,  z£Q), 

G(z,  z)  =  -^-^+(l-4v+8v2)[R,{2z)+ilog(z/4a)] 

IT  C  TT 

+  ~2  [(l~4v)azR-(2z)+z2R4(2z)](0<z$t),  G(0,  04)=  -  4v(1+4v)  t 
a  w  / 


}  (2.23) 
f  cont. 


Equations  (2.  19).  (2.  22),  (2.  23),  together  with  the  elliptic- 
integral  representation  of  R^  (n=2,3,4)  contained  in  (2.  16),  supp’y  —  in 
a  form  suited  to  numerical  purposes  —  the  kernels  8f(z,  Q/&C  and 
9g(z,  C)/8C  involved  in  (1.  27),  (1.  28)  for  the  case  of  a  bar  of  circular 
cross-section.  Further,  the  values  f(z,  0),  g(z,  0)  (O^z^-t),  which  are 
also  needed  in  connection  with  (1.  27),  (1.  28),  are  computable  from 
(2.  11),  (2.  16)  for  0<z£<-,  while 

f(0+,0)  =  l,  g(0+,0)=(l-2v)2,  (2.  24) 1 

as  is  easily  verified  on  the  basis  of  (2.  11),  (2.  18).  The  Fredholm 
integral  equation  (1.27),  which  determines  the  fictitious  bar-force 

2 

P#(z)  (OSzrJ-t,),  may  now  be  solved  numerically  for  any  admissible  choice 
of  the  elastic  constants  V,  13,  and  v.  To  accomplish  this  task  the  range 
of  integration  [0,4.]  in  (1.  27)  was  partitioned  uniformly.  The  contribu¬ 
tion  to  the  improper  integral  arising  from  the  (integrable)  logarithmic 

3 

singularity  of  the  kernel  8g(z,  Q/dQ  at  z=Q  was  then  evaluated  in  closed 
elementary  form  (in  terms  of  the  values  of  p+  at  the  mesh-points  of  the 
partition)  upon  replacing  p+  by  a  continuous,  piecewise  linear  function. 

On  the  other  hand,  the  ordinary  trapezoidal  rule  was  used  to  compute 

^  Note  the  agreement  between  the  first  of  (2.  24)  and  the  last  of  (1.  30). 

2^  Recall  from  (1.  4)  that  r^rr^-r^O  and  see  the  footnote  attached  to  (1.  14). 

3) 


See  (2.  19). 
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the  contribution  to  the  integral  under  consideration  stemming  from  the 
bounded  portion  G  of  the  kernel,  which  is  given  by  (2.  23).  In  this 
manner  (1.  27)  was  reduced  to  a  system  of  linear  algebraic  equations, 
whose  solution  was  obtained  on  an  electronic  computer.  Once  p^z) 
(0£z^t)  had  been  so  determined,  the  desired  actual  bar-force  p(z) 
(0£z£t)  was  found  from  (1.  28)  by  means  of  a  numerical  evaluation  of 
the  improper  integral  in  (1.  28)  that  is  strictly  analogous  to  the  numer¬ 
ical  integration  scheme  described  above. 

Before  proceeding  to  the  results  thus  established,  we  investigate 
the  asymptotic  behavior  of  p(z)  at  the  endpoints  z=0  and  z-l  of  the  em¬ 
bedded  bar-segment.  Substituting  for  o  from  (1.  24)  into  (1.  13)  and 
differentiating  the  resulting  identity  with  respect  to  z,  one  draws 


p<*>4{[p0-p«<o>}!§|(z  0)+p*(t>!i 


(z,t) 
l 

-Ip*(c)fl£(z’£)dc}  (°<z< l>- 


\ 

>  (2.25) 

/ 


Moreover,  for  the  bar  of  circular  cross-section,  8f(z,()/8z  admits  the 
representation  (2.20).  Next,  we  suppose  that  for  some  fixed  a>l,  p^ 
is  absolutely  integrable  on  [0,-t]  and  invoke  the  Holder  inequality*  to 


justify  the  estimate 


0  0  0 


tt>l,  -+^  =  1,  (0$Z£‘t). 


TT 


See,  for  example,  Beckenbach  and  Bellman  [13],  p.  21. 


\ 


>  (2.26)2 


) 


21 

'  Note  that  the  second  integral  in  the  right-hand  member  of  (2.26)  exists 
by  virtue  of  (2.20),  for  every  (3>0.  Further,  it  is  not  difficult  to  show 
that  this  integral  is  bounded  for  0£z£t. 
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Consequently,  the  integral  in  the  right-hand  member  of  (2.25)  is 
bounded  on  (0,  l)  under  the  foregoing  rather  weak  assumption  concern¬ 
ing  p^.  In  these  circumstances  it  follows  from  (2.  25)  and  (2.  20)  that 

p(z)=0(l)  as  z-»0,  'l 

(l-2v)P)jl(-t)  >  (2.27) 

P(z)=  2n'(r-v) a  -lo8^H-z)/a]-»-0(n  as  z-t  J 

Next,  it  is  essential  to  examine  the  manner  in  which  the  various 
physical  parameters  enter  the  solution  of  the  problem  under  consider¬ 
ation.  To  focus  attention  on  this  issue,  we  now  write 

P(z)=p(z;po,  a,  t,  t/,  r|,  v)  (Oizzl),  (2.28) 

where,  it  will  be  recalled,  the  parametric  arguments  pQ,  a,  l,  r r|  and 
v,  in  this  order,  stand  for  the  applied  scalar  bar-load,  the  radius  of 
the  bar,  the  length  of  the  embedded  bar-segment,  Young's  modulus  of 
the  bar.  Young's  modulus  of  the  surrounding  medium,  and  Poisson's 
ratio  of  the  latter.  It  is  clear  from  the  structure  of  (1.27),  (1.28),  in 
view  of  (1.4),  (2.  11),  (2.  13),  and  (2.  10),  that  the  solution  may  be  cast 
into  the  dimensionless  form 

p(z)/po=(r(z/a;  Va,  r\lr\,  v)  (Osz/a^t/a),  (2.29) 

so  that  the  dimensionless  (normalized)  bar-force  p(z)/pQ  depends  exclu¬ 
sively  upon  the  dimensionless  depth-coordinate  z/a,  the  length-ratio  -t/a, 
the  stiffness -ratio  r\/r\,  and  the  Poisson- ratio  v.  The  function  ♦  in 
(2.  29)  is  independent  of  l/a  and  v  for  T],/'n=l,  i.  e.  ,  if  the  bar  and  the 
embedding  medium  have  the  same  modulus  of  elasticity.  For  this 
degenerate  limiting  case,  (1.  32),  the  first  of  (2.  11),  and  (1.25)  furnish 

p(z)=2poCR1(z)+^R2(z)]  (0<zst),  p(0)=po,  (2.30) 


with  Rj,  R2  given  by  (2.  15),  (2.  16). 
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We  turn  now  to  the  discussion  of  illustrative  numerical  results 
for  the  decay  of  the  actual  bar-force  p,  which  are  depicted  graphically 
in  Figure  3  to  Figure  8.  The  load-diffusion  curves  displayed  in  these 
diagrams  show  the  normalized  bar-force  p(z)/pQ  as  a  function  of  the 
dimensionless  depth-coordinate  z/a  for  O^z/as-t/a  and  various  choices 
of  the  dimensionless  physical  parameters  Va,  V/t),  v>  Each  of  the  six 
figures  under  discussion  refers  to  a  single  Poisson's  ratio  v  of  the 
embedding  medium,  to  two  values  of  the  length-ratio  li a,  and  to  four 
distinct  choices  of  the  stiffness-ratio  T\ /r\.  The  thirty-six  combinations 
of  parameter- values  covered  by  the  graphs  presented  here  are  listed 
in  the  following  table. 

Figure  3:  v=0;  £/a=  5,10;  V/t|=1, 2,  4,  8. 

Figure  4:  v=0;  £/a=10,20;  T|7t)  =  1, 2,  4,  8. 

Figure  5:  v=^;  Va=  5,  10;  V/n=l,  2,  4,  8. 

Figure  6;  v=^;  -t/a=10,20;  n7n=l,  2,  4,  8. 

Figure  7:  v^;  Va=  5,  10;  T//n=l,  2,  4,  8. 

Figure  8;  v=j;  Va=10,20;  r|#/r)=l,  2,  4,  8. 

As  is  apparent  from  this  table,  the  load-diffusion  curves  appropriate 
to  Va=10  in  Figures  3,  5,  7  are,  for  comparison  purposes,  repeated 
in  Figures  4,  6,  8. 

All  of  the  curves  under  consideration  represent  functions  p(zJ/pQ 
that  are  steadily  decreasing  throughout  their  interval  of  definition 
0£z/a£{/a.  Further,  in  each  instance, 
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P<0Vpo  =  1.  p(t)/po>0,  (2.31)1 

the  first  of  (2.  31),  according  to  (1.  31),  being  valid  for  a  bar  of 
arbitrary  cross-section.  Therefore,  no  part  of  the  applied  load  is 
communicated  to  the  embedding  medium  through  bond-forces  concen¬ 
trated  at  z  =  0,  whereas  a  portion  of  the  load  is  transmitted  to  the 
surrounding  solid  through  bond-forces  concentrated  at  the  terminal 
cross-section  of  the  bar.  The  graphs  indicate  that  this  load-portion 
increases  with  the  stiffness -ratio  and  is  a  decreasing  function  of  the 
length-ratio;  in  the  examples  to  which  Figures  4,  6,  8  refer,  the  frac¬ 
tion  of  the  load  transferred  at  the  embedded  end  of  the  bar  is  less 
than  10%  for  -t/a=l0  and  less  than  5%  for  Va=20. 

The  absence  of  any  concentrated  load-transfer  at  z=0,  predicted 
by  the  current  solution,  is  interesting  in  view  of  the  corresponding  re¬ 
sult  encountered  earlier  in  an  analogous  treatment  of  Reissner's  plane 

2 

load-transfer  problem  ;  there,  a  substantial  part  of  the  load  was  found 
to  be  transmitted  to  the  sheet  at  the  junction  of  sheet  and  stringer.  On 
the  other  hand,  as  was  shown  in  [5],  singular  load-transfer  at  either 
end  of  the  attached  stringer-segment  is  precluded  within  Reissner's  [2] 
original  formulation  of  the  plane  problem. 

For  all  illustrative  examples  considered  here,  the  initial  slope 
(at  z=0)  of  the  load-diffusion  curve  is  finite  in  agreement  with  the  first 
of  (2.27);  the  infinity  of  the  end-slope  at  z=t,  predicted  by  the  second 

^  The  analytical  determination  of  p(t),  on  the  basis  of  (1.  27),  (1.  28), 
appears  to  offer  serious  difficulties  even  for  the  circular  bar.  In 
the  degenerate  case  p(t)  for  the  circular  bar  is  immediate 

from  (2.  30),  (2.  15),  (2.  16). 


2) 


See  Section  3  of  [5]. 
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of  (2.27)*  if  v^l/2,  is  not  discernible  on  the  scale  of  Figure  3  to 
Figure  6. 

As  is  intuitively  plausible,  each  of  the  diagrams  reflects  a  less 
rapid  diffusion  of  the  load  with  rising  values  of  the  stiffness -ratio. 

Even  for  a  relatively  stiff  bar  of  sufficiently  large  length- ratio,  however, 
most  of  the  load-transfer  takes  place  in  the  vicinity  of  the  surface  of  the 
embedding  medium  and  —  in  this  range  —  is  quite  insensitive  to  changes 
of  the  length- ratio.  To  illustrate  this  observation,  we  refer  to 
Figure  6  (v=l/4)  and  regard  the  bar- radius  as  fixed.  An  inspection  of 
the  curves  in  this  figure  reveals  that  for  ,n//r|<8  and  -t=20a  over  half  of 
the  applied  load  is  transferred  to  the  matrix  by  the  bar- segment  in  the 
range  0<z£4a.  Further,  the  dashed  curves  appropriate  to  1=1  Oa  are, 
in  this  range,  practically  indistinguishable  from  the  corresponding 
solid  curves^,  which  pertain  to  <t=20a.  Accordingly,  here  very  little 
benefit  is  derived  from  a  doubling  of  the  embedded  bar-length.  This 
finding,  which  is  also  borne  out  by  the  results  given  in  Figure  4  and 
Figure  8,  would  appear  to  be  significant  for  design  purposes. 

In  conclusion  we  emphasize  onre  more  that  the  analysis  carried 
out  in  this  section  is  bound  to  be  inadequate  unless  the  length  of  the 
embedded  segment  of  the  bar  is  sufficiently  large  compared  to  the  bar- 
diameter. 


P#W  was  *ound  not  to  vanish  in  the  examples  dealt  with  numerically. 

The  dashed  and  the  solid  curve  for  r|7n=l  arc  identical  in  all  figures 
since  here  p  is  independent  of  the  length- ratio,  as  well  as  of  v. 

See  (2.  30)  and  the  remarks  preceding  (2.  30). 
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Appendix;  Proof  of  the  limit- relations  (1.  23). 

We  indicate  here  a  proof  of  the  limit- relations  (1.23)  pertaining 
to  the  functions  (n=l,2,3)  defined  by  (1.  19),  (1.  20).  In  order  not  to 
detract  from  the  essence  of  the  proof,  we  shall  confine  our  attention  to 
a  cross-sectional  domain  II,  the  boundary  811  of  which  is  a  simple 
closed  curve  with  continuously  varying  curvature *.  Accordingly, 
denoting  by  s  the  arc -length  of  dll  (measured  counter-clockwise  from 
an  arbitrarily  chosen  fixed  point  of  the  closed  curve  at  hand)  we  assume 
that  8 II  admits  the  parameterization 

811:  XpXjts),  x2=X2(s)  (0sb<so),  (1) 

where  Xj  and  X2  are  functions  twice  continuously  differentiable  on 
[0.8  ].  while 

y0)=xp(8o).  yo)=ip(so),  yo)=yso)  0=1,2).  (2) 

sq  being  the  total  length  of  811.  The  curvature  H  of  8ll  is  given  by 

h(s)=X1(s)X2(s)-12(s)X1(s)  (3) 

and,  for  convenience,  we  set 

m=  max  |*(s)|,  t  =^-.  (4) 

O^sSs  ° 

o 

Next,  let  II'  be  the  closed  band-shaped  sub-region  of  II  consisting 
of  all  points  in  H  whose  perpendicular  distance  from  811  fails  to  exceed 

tjj,  ©•  , 

n/  =  {(x1,x2)|(x1,x2)€ff,  [x1-X1(s)]2+[x2-X2(s)]25t^  (0ss5so)}.  (5) 

Further,  for  every  point  (x^,x2)  in  II/,  let  s(x^,x2)  be  the  arc-length 

)  The  following  argument  is  easily  generalized  to  accommodate  a 
boundary  8ll  that  exhibits  corners,  but  is  composed  of  a  finite  number 
of  arcs  with  continuous  curvature. 
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associated  with  that  orthogonal  projection  of  (x^,  x^)  onto  dll  which  has 
the  smallest  distance  from  (x^x^);  also,  let  t(xj,x2)  be  the  distance 
of  this  orthogonal  projection  from  the  point  (x^.x^).  Evidently, 
s(Xj,x2)  andt(xj,x2)  define  an  orthogonal  curvilinear  coordinate-net 
on  II'.  The  underlying  coordinate  transformation  is  characterized  by 


x1(s,t)=X1(s)-tX2(s), 


(6) 


x2(®»  (0^8<so#  0^t^tQ)f  j 

and  the  mapping  (6)  is  one-to-one  on  II*;  indeed,  its  Jacobian  obeys 
0x.  0x?  0x.  0x_ 

J(8*t>  =  S^5F''BF"57_  =  1"tH<8>>0  (0^8<so,0etsto),  (7) 

the  inequality  being  a  consequence  of  (4).  Finally,  for  every  (x^.x^) 
in  n,  we  adopt  the  notation 


S(x1.x2)  =  {(?1.  ?2)I(?1-x1)2+(?2-x2)2sM2(x1.x2)}. 

U(xl,x2)=minV(?1-x1)2+(?2-x2)2,  (5j,  ?2)G0lI, 


>  (8) 


so  that  E(Xj,x2)  stands  for  the  interior  of  the  largest  closed  circular 
disk  contained  in  I!  and  centered  at  (X|,x2). 

Having  disposed  of  these  geometric  preliminaries,  we  recall 
from  (1.  20),  (1.  19),  and  (1.  3)  that 


where 


wl<*)4  Vl(x1,x2,z)dA  (0<  |z|<oo). 


(9) 


1 


TT 


Note  that  we  now  write  (x^.x^x^)  in  place  of  the  single  vectorial 
argument  x. 
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V1(x1.x2,z)  =  ^  J(p2+z2f3/2dA,  ,  P2=(x1-§1)2+(x2-?2)2  (10)1 


for  all  (xj,x2)€I1  and  all  z^O.  From  (9)  follows  the  estimate 

lwl(z)-^h<z)|*M(z)+N(z)  (0<|z|<co),  (11) 

in  which 

M(z)  =  -^  J  |vi(x1,x2,  z)-^fh(z)|dA, 
n-n' 

/  (12) 

=  J  |vi(x1,x2,z)-^h(z)|dA  (0<|z|<oo). 
n' 


whereas  h  is  the  step-function  defined  in  (1.  25).  In  order  to  establish 
the  first  two  of  (1.  23)  it  suffices  to  show  that  M(z)  and  N(z)  tend  to 
zero  as  z-*0. 

Now  fix  (x^,x2)  in  II,  write  E  and  |i  in  place  of  E(xj,x2)  and 
U(xj,x2),  and  infer  from  (10)  that 


U 


.  .  2trz  f  pdp  ,  z  f  ._2,  2.-3/2ja 

V1(x1,x2,z)  =  ir  J  2  2  5/Z  +  j  J  (P  +  z  )  dA 


0  (p  +z  )  n-E 


for  every  z£0.  However, 


z  I  S  =h<z>-  J==~ 
o<“  >  /U2+I2 


CD 


l  (p2«2,-3/2dA?<2nJ  PdP 


2ir 


n-E  u 

Combining  (13)  and  (14),  one  has 


777707? -TFT 

~  u  (p  +z  >  /u2+z2 


(z*0). 


(13) 


A 

/  (14) 


/ 


TT 

2) 


The  subscript  5  attached  to  the  element  of  area  is  to  indicate  that 
(?1»  ?2)  are  the~ variables  of  integration. 

See  (8)  for  the  definition  of  E(x1#x2)  and  u(xj,x2). 
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lvl<xl*x2*  z)"'Xh^z^< 


4rr  z 


(0<|z|<cd). 


(15) 


a/i2+z*' 

But,  according  to  (5)  and  (8), 

0<tQ<|i(Xj,  x^)  for  all  (xj,  x^)^!! -II/.  (16) 

Consequently,  if  A*  stands  for  the  area  of  II#,  (15)  and  the  first  of  (12) 
furnish 


M(z)<1?(4-a')I*Lo(1) 


A2/2^2' 
o 


as  z-»0. 


(17) 


On  the  other  hand,  the  second  of  (12),  together  with  (15)  and  the 
second  of  (8),  yield  the  estimate 


N(z)<±^Ll(z),  I(Z)=J 


dA 


'  fJ, 

/t  (Xj.X. 


)+z 


21 


(0<|z|<oo) 


(18) 


since*  u(Xj,  x^tfxj,  x^)  for  ail  (x^.x^)  in  II#.  Changing  the  variables  of 
integration  in  (18)  from  (x^.x^)  to  the  curvilinear  coordinates  (s,t)  intro¬ 
duced  in  (6)  and  bearing  (7)  in  mind,  we  arrive  at 


t  8 
Q  o 


«*>=  J  I  dsdt  (0<|z|<oo). 

0  0 


(19) 


From  (19)  and  (4),  in  turn,  follows  the  estimate 
t 


0<I(z)s 


2S  3  s 

°  r  dt-=  °[log(t  +  /t2+zZ’)-log|z|]  (0<|zj<oo).  (20) 

fa  o  o 

z 


J  — 


Combining  (18)  and  (20),  we  reach 


N(z)<4TrLZl  I(z)=o(l)  as  z-*0. 
AT 


(21) 


TT 


Recall  the  geometric  significance  of  the  curvilinear  coordinate 
t(Xj,X£)  associated  with  the  mapping  (6). 
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Finally,  (11).  (17),  and  (21)  assure  that 

lim[w1(z)-^h(z)]=0,  (22) 

z-*0 

which  confirms  the  first  two  of  (1.  23).  The  remaining  limit- relation 
in  (1.  23)  may  be  established  by  similar  means. 
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FIGURE  I.  GEOMETRY  OF  BAR  AND  EMBEDDING  MEDIUM 
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FIGURE  2.  BAR,  EXTENDED  EMBEDDING  MEDIUM,  FICTITIOUS  REINFORCEMENT 
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FIGURE  3. 
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FIGURE  4.  CIRCULAR  BAR  v  =  0,l/a 


FIGURE  5.  CIRCULAR  BAR’  v  *  1/4,  1/a  *  5,  //o  » 10 
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